MECHANICS OF COMPLEX BODIES: COMMENTARY ON THE 
UNIFIED MODELLING OF MATERIAL SUBSTRUCTURES 
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Abstract. Basic issues of the general model-building framework of the me- 
chanics of complex bodies are discussed. Attention is focused on the repre- 
sentation of the material elements, the conditions for the existence of ground 
states in conservative setting and the interpretation of the nature of the various 
balance laws occurring. 
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1. Introduction 

Materials in which changes in the molecular or crystalline texture at various mi- 
croscopic scales (substructure) influence the macroscopic behavior through peculiar 
interactions are commonly available. Liquid crystals, ferroelectrics, quasicrystals, 
polymeric fluids are paradigmatic examples. The attribute complex is assigned to 
bodies made of these materials in order to underline that significant substructural 
effects must be accounted for. 

In complex bodies the prototype material element is a system. Often it is a 
perfectly identifiable Lagrangian system, like in nematic liquid crystals in which 
the characteristic stick molecules embedded in a soft matrix may be extracted 
and isolated from the rest. Sometimes it does not and the substructure is in a 
certain sense virtual, like in microcracked bodies: microcracks do not exist per se, 
rather they are determined only by the surrounding matter. On the other hand, 
the substructure can be procedural in the sense that it is a consequence of local 
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rearrangements of the matter due to phase transitions from one energetic well to 
another as in martensite-austenite mixtures. 

Notwithstanding the variety of phenomena displayed by complex bodies and clas- 
sified in condensed matter physics, there exists an abstract model-building frame- 
work for the mechanics of complex bodies. It unifies in a unique format existing 
models of special classes of complex bodies and is a flexible tool for analyzing new 
materials. 

Basic aspects of such a model-building framework are discussed here with peda- 
gogical purposes and constructive criticism. Some appropriate references are scat- 
tered throughout the subsequent sections. 

2. Representation of the morphology of material elements 

In its primitive meaning, a body can be regarded as an abstract set B collect- 
ing material elements, each one being the smallest piece of matter characterizing 
the material composing the body. The basic issue is the 'representation' of such a 
set. In the standard format of continuum mechanics the geometrical representation 
adopted is the minimal one: each material element is mapped onto a place that it 
occupies in the ambient space R d . However, in the real world the material elements 
are groups of entangled molecules, simple or complex pieces of crystalline struc- 
tures, stick molecules dispersed in a ground fluid etc., and substructural changes 
may generate actions that it is hard to identify only with perturbations of the 
standard stress. In all these cases the standard representation of bodies is too min- 
imalist. The material element should be considered in essence a system rather than 
a windowless box so that, in the representation of the body, a map attributing to 
each material element a morphological descriptor v of its (inner) substructure has 
to be defined. To construct the essential structures of the mechanics of complex 
bodies, at least at the level of first principles, it is not necessary to render precise 
the nature of the morphological descriptor except assuming that it is an element of 
a finite-dimensional differentiable manifold M. with minimal geometrical properties 
(each property of M. has in fact physical meaning, see [H]). Once a reference place 
B for the entire body is selected in M. d , any other actual place B a (with the index a 
meaning actual) is considered to be achieved in an isomorphic copy M. d of R d itself 
by means of a transplacement map 

B 3 x ^ y (x) £ B a 

which is assumed to be one-to-one, differentiable and orientation preserving, with 
spatial derivative 

F := Du (x) e Horn (T X B, T u{x) B a ) . 

The geometry of the inner structure of the material elements is described (at least 
at a coarse grained level) by a morphological descriptor map 

B 3 x i — > v(x) £ M 

which is assumed differentiable with spatial derivative 

N := Dv (x) £ Horn (T X B, T u(x) M) . 

Motions are then time-parametrized families of transplacements and morphological 
descriptors, namely for x £ B and (G [0,(] cl the time 
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(x, t) i — > v := v (x, t) E M, 

both fields assumed to be piecewisc twice diffcrcntiablc in time so that 

d , , . d . . 

y-=^y{x,t) and v.= —v[x,t) 

represent the macroscopic velocity and the rate of change of the substructure re- 
spectively. 

3. Standard and substructural actions 

3.1. Observers. The description of material substructures constrains to give a de- 
tailed look to the essence of the notion of observer. I have discussed the point 
repeatedly (see e.g. [12] and references therein). The essential aspects of the dis- 
cussion are summarized below (with additional remarks) as a preamble to the use 
of a weaker invariance requirement under changes in observers presented later. The 
aim is to underline how the geometrical features of the ambient space R d and the 
manifold of substructural shapes M. influence the structure of the integral balances 
of the interactions of macroscopic and substructural nature. 

My point of view is that an observer should be considered as a representation of 
all geometrical environments necessary to describe the morphology of a given body 
and its motion. To this aim, in standard continuum mechanics one needs to select 
the reference place B, the ambient space, say R d , and the interval of time. They are 
all the geometrical environments needed. Different is the stage when substructural 
complexity arises and its changes influence the gross behavior. The manifold of 
substructural shapes comes into play in the way described above and one needs to 
represent it, specifically, one needs to select atlantes over M. 

Once M is accounted for, the definition of changes in observers should involve the 
representations of R d , £>, [0, t] , M.. In particular, the attention is focused here on 
changes in observers leaving invariant B and the interval of time. The last require- 
ment defines synchronous changes in observers. Really, one could consider affinc 
time rescaling: non essential consequences accrue so that they are not considered 
here. 

• Changes in the ambient space are given by elements of Diff(R d ,R d ), the 
group of diffcomorphisms of R d onto itself. Really one takes smooth curves 

si — ► f s £ Diff(R d , R d ), seR+, 

with fo = id, where id means identity. The parameter s can be identified 
with the time so that the curve s i — > f s can be interpreted in the common 
way as the motion by which two observers differ as the time flows. In 
particular the vector field y i — ► v (y) :— 4^ | s= o can be considered as 
the (virtual) velocity of an observer moving with respect to another one, a 
velocity pulled back in the frame of the first observer. 

• The material substructures are placed in the ambient space and the man- 
ifold M. collects only the elements of a concise description of the charac- 
teristic features of their geometry in space. Thus changes of frames in the 
ambient space alter in principle the geometry of the substructures and their 
consequent representation over M.. Disconnection between changes in the 
manifold of substructural shapes and changes of frames in space is admissi- 
ble only when v represents only a generic property of the substructures not 
associated with their geometry in space. Changes in the choice of atlantes 
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over M. are governed by elements of the Lie group of diffeomorphisms of 
M. onto itself, namely 

G := {g : A4 — > M. \ g a diffeomorphism} . 

The link with changes in the ambient space are then assured by assuming 
the existence of a homeomorphism 

h : Diff{R d ,R d ) -> G, h(id)=id G , 

where id,Q is the identity over G, so that a curve 

s i — ► ^ 3 := h s (v) , 

with h s = h(f s ). By indicating by £ the element of the Lie algebra g of G 
given by the derivative ^ | s= o = ^jj^ |s=o , its value over a given v is in- 
dicated by $, M (u) := ^p- | s= o . In particular, if the curve s i — ► i s is selected 
over the special orthogonal group SO (d), a subgroup of Dif f(M. d ,R d ), for 
Ag an element of the Lie algebra so (d), q € M. d , it is possible (and also conve- 
nient) to write ^ M (y) as the product A (v) q with A (v) £ Hom(R d , T V M). 
Isometric semi- classical changes in observers are the ones that, by leaving in- 
variant B and [0, t] , are characterized by the choice of the subgroup of Diff(M. d , M. d ) 
coinciding with the semi- direct product R d tx SO(d). 

Rotational semi- classical changes in observers are the ones in which the subgroup 
of Dif f (Mr, R d ) selected as ambient of s i — ► f s is just SO (d). In this case, by 
indicating by y* and v* the rates evaluated after the change in observer (they are 
the pull-back in the frame of the first observer of the rates evaluated by the second 
observer), one gets 

(3.1) y* = y + qA(y-y ), 

where yo is an arbitrarily fixed centre of rotation in the ambient space, and 

(3.2) v* = v + A [v) q. 

Remind that q depends only on the parameter s that is identified here with the 
time, so that q = q (t). 

3.2. Augmented external power and SO (d) invariance. A part of B is any 
subset b of B itself with non-vanishing volume measure and the same geometric 
'regularity' of B (that is the same topological properties). All parts of B form an 
algebra *p (B) with respect to the operations of meet and join (see [5])- 

Let Vel be the space of all rate fields (x, t) i — ► y (x, t) and (ar, t) i — ► v (x, t) over 
the tube B x [0,i], rates calculated along possible motions (x,t) i — > y(x,t) and 
(x, t) i — ► v (x, t). 

A power along motions (y, v) is a functional 

V : «P (B) x Vel -> R 

which is additive over disjoint parts and linear in the rates. 

The explicit representation of a given P requires classification of the interactions 
occurring in a body. The interest here is on the expression of the power of all exter- 
nal actions over a generic part b and along a motion (y, v), a functional indicated 
from now on by (y, v). 

Once b is selected arbitrarily, interactions with the rest of the power and the ex- 
ternal environment are classified in two main subclasses: (1) standard actions power 
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conjugated with the macroscopic deformation, (2) substructural actions associated 
with the rate of change of the substructure inside the material elements. Each 
subclass is further subdivided into bulk and contact actions which admit densities 
with respect to volume and surface measures dx and d7i d ~ 1 , respectively. 

The natural expression of the external power satisfying previous assumptions is 
then given by 

Vt xi {y,P):= [ (b-y + P-v) dx+ [ (t-u + T-v) dH d -\ 

Jb Jdb 

a power written in referential form. The term 'augmented' in the title of this section 
underlines the presence of the power densities of the substructural actions. At any 
x in b one gets 

b := b (x) e T; (x) B a ~ R d , 0:=0 (x) e K (x) M, 

while, for x G db, 

t := t (x) £ T* y(x) B a ~ R d , r r (x) e T* u{x) M. 

Cauchy theorem indicates that the standard traction t can be expressed in La- 
grangian representation by means of the first Piola-Kirchhoff stress tensor P, so 
that t = Pn at any x G db where n is defined. P (x) belongs to Hom(R d ,T* (x) B a ) 
and also it is natural to consider it as the density of a form over B. The argument of 
the proof is the standard Cauchy's one based on the tetrahedron, or subtle refine- 
ments of it (see the results in |20|). In analogous way one may presume that r = Sn 
with S (x) 6 Hom(M. d , T*^ M). Of course, even the microstress S can be consid- 
ered as the density of a form over B. The standard tetrahedron argument does 
not apply here because the field x i — > r (x) takes values on the whole cotangent 
bundle T* M. = U ve MT*M so that its total over a generic side of the tetrahedron 
is not defined, as pointed out more in general later. To prove the existence of S by 
using the tetrahedron argument, it is necessary to embed M. in a linear space (the 
relevant proof is in [5]). The embedding always exists since M. is assumed here to 
be finite-dimensional (Whitney theorem) and can be selected to be also isometric 
(Nash theorem), this choice having the advantage to preserve the quadratic part of 
the independent kinetic energy that can be sometimes attributed to the material 
substructure (see relevant comments in |13j). However, the theorems indicating the 
availability of the embedding of M. in a linear space do not assure the uniqueness 
of the embedding itself. More precisely, the embedding is neither unique nor rigid 
(A4 can be at the end folded in various manners in the process). Additionally, the 
dimension of the target linear space depends on the regularity of the embedding 
(Nash theorem). In all cases, if one finds convenient embedding A 7 ! in a linear space 
for technical purposes, for example for constructing appropriate finite elements for 
some special model of complex bodies, the choice of the embedding becomes strictly 
a matter of modelling. Of course, such a peculiarity disappears when the complex 
material under examination admits a manifold of substructural shapes which is 
coincident with a linear space, as in the case of micromorphic or, more generally, 
affine bodies (see [18]). 

One might naively claim that M. is always a linear space so that the scheme of 
affine bodies (the one discussed for example in [7], [9]) is sufficient for analyzing 
the material complexity at substructural level. 
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Solids with distributed magnetic spins in conditions of magnetic saturation (M. 
coincides with S 2 ) and the superfluid helium A He (where M. = S 1 C C) are ele- 
mentary counterexamples to the claim. With the aim of unifying the treatment of 
as many special cases of physical interest as possible, it is then necessary to con- 
sider A4 as an abstract manifold prescribing only the minimal geometrical prop- 
erties necessary to build up the essential objects which are useful for constructing 
the mechanics of complex bodies. In this case the existence of the microstress 
S (x) G Hom(M. d ,T*^A4) can be assumed a priori (an intrinsic representation of 
S in terms of measures is presented in |19j where the primary object considered is 
the inner power, instead of the external one). 

Without investigating further on the question, I assume here that t depends 
linearly on the normal at x to db so that the natural expression of the external 
power of all actions over the generic part b along y and v is then 

Vt xt {y,v):= [ (b-y + p-u) dx + [ (Pn ■ y + Sn ■ v) dH^ 1 . 

Jb Jdb 

A crucial axiom is the requirement that the power V^ xt is invariant under iso- 
metric changes in observers, that is under the action of the semi-direct product 
M. d tx SO (d) over the ambient space and the action of elements of G over M. in- 
duced by the homomorphism h introduced above (see also additional remarks in 
[llj . [13] , [15]). Here, a weaker axiom is used, namely invariance of the power un- 
der the sole action of SO (d) in space and the corresponding action of G over M. 
through h. It is thus required that observers differing only by a proper rotation 
evaluate the same power. 



Axiom 1. (SO (d) invariance) At mechanical equilibrium the external power of all 
actions on any part of B is invariant under rotational semi- classical changes in 
observers, namely 

(3.3) vr t {y\o*) = vr t {y,u) 

for any choice of the rotational velocity q (t) £ R d and any b £ (B) . 



Theorem 1. (i) If for any b the vector fields x i— > Pn and x A*Sn are defined 
over db and are integrable there, the integral balances of actions on b hold: 

(3.4) [ bdx+ [ Pn dH^ 1 = 0, 

Jb Jdb 

(3.5) f {(y-yo) Ab + A*(3) dx + [ ((y - yo) A Pn + A*Sn) dTL^ 1 = 0. 

Jb Jdb 

(ii) Moreover, if the tensor fields xi->P and x i— ► S are of class C 1 (Bo) fl C°(23 ) 
then 

(3.6) DivP + 6 = 

and there exist a covector field x i— » z G T v i x \M such that 

(3.7) skw{PF*) = e(A*z + (DA*)S) 
and 

(3.8) DivS -z + P = 0, 
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with z = z\ + Z2, Z2 G KerA*. (Hi) If the rate fields (x,t) i — ► y(x,t) G M. d and 
(x,t) i — > i>(x,t) G T v mM. are differ entiable in space, the local balances imply 

(3.9) Vt xt {y,v) = vt t (y,y) 

where 

rl nt (y, z>) := I \p ■ F + z ■ v + S ■ N) dx. 

Jb 

Above e is Ricci's alternating index. V™ 1 (y, v) is called an inner (or internal) 
power. 

Proof. The immediate consequence of the axiom of SO (d) invariance is the integral 
balance of moments (|3.5j) obtained by using ()3.1|) and (|3.2|) in (|3.3p . In (13. ip the 
point yo is arbitrary. As a consequence, by taking an arbitrary vector w G M. d , one 
can substitute yo with yo + w in (|3.5p . Such a substitution corresponds to a simple 
shift of the centre of the rotation of one observer with respect to the other. By 
subtracting (|3.5p from its counterpart calculated at yo + w, one then gets 

w([bdx+[ PndH d - 1 ) = 0, 

Jb Jdb 

which corresponds to (|3.4[) as a consequence of the arbitrariness of w. Note that 
the substitution yo — > yo + w is possible due to the linear structure of M. d , a 
structure that is in general not available over the manifold of substructural shapes 
M. . Under the regularity hypotheses above, the local balance (|3.6|) follows as usual 
by exploiting Gauss theorem and the arbitrariness of b. The same localization 
procedure applied to the integral balance (|3.5|) and the validity of (|3.6p imply the 
local balance 

ePF* - (DA*)S = A* (DivS + (3) . 

Since A* {v) G Hom(T*M.,M. ), two information are available from this equation: 
(1) At each v G M the difference ePF* - (DA*) S is the image in R d of a covector 
in T*M, let say z. (2) Such a covector is just equal to DivS + [3. The equation 
(|3.9p follows by direct calculation under the validity of the pointwise balances Q3.6P 
and (3J§. □ 

Of course, the balance equations above include the dynamic case because the 
bulk actions can be decomposed additively in their inertial and non-inertial parts, 
the latter being identified by requiring that their power balances the rate of change 
of the kinetic energy. In this procedure it is assumed that the energy is the sum of 
macroscopic and substructural contributions. For the sake of conciseness the topic 
is not developed here (see [13] for the details). 

Theorem 1 is the same that can be obtained by imposing as an axiom the in- 
variance of the external power under isometric semi-classical changes in observers 
mentioned earlier (see jllj). Such an equivalence underlines that the integral bal- 
ance of standard forces is a peculiar consequence of the 'rigid' structure of R d (its 
linear structure) and is associated with one of the Killing fields of the metric in 
the ambient space (see [H] for the analogous observation in the case of simple 
bodies). The same property is not available over Ai straight away. In fact, it has 
been assumed here that the manifold of substructural shapes is abstract so that 
it does not coincide with a linear space in general. For this reason the totals of 
the substructural actions are not defined a priori, as mentioned above in discussing 
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the possible path toward the proof of the existence of the microstress S. Consider, 
for example, the field x i — ► (3 (x) G T*,sM that is : B — > T*M. The target 
space T* M. is not a linear space so that the integral of (3 on any part b of B is 
not defined unless M. itself is a linear space. Analogous remarks hold for the fields 
x i — ► z (x) G T*,sM. and x i — > (<Sn) (x) G T^^AW. As a consequence, not only 
an integral balance of substructural actions does not follow from the requirement 
of invariance of the power under changes in observers but it is even not defined. 
Moreover, the fact that j3 and S appear only in the balance of moments does not 
means that they represent (micro) couples because in the integral balance of mo- 
ments they are multiplied by the the formal adjoint of A which maps at each v 
elements of T V M. onto elements of R d . 

Various alternative paths can be followed to get pointwise balances of actions. 
The comments below apply to them. 

(1) One could postulate the integral balances of standard and substructural 
actions as first principles. However, as pointed out above, such a point of 
view can be adopted only in the (very) special case in which AW is a linear 
space. When it is not the case, the balance of substructural actions cannot 
be postulated, because the integrals appearing are not defined. 

(2) One could adopt the virtual power procedure proposed in [7| for affinc 
bodies, by postulating in fact the weak form of the balance equations. In 
this case one must postulate not only the expression of the external power 
but also the internal power, the power of the inner actions. In this way 
one should postulate the existence of the inner self-action z. In contrast, 
in Theorem 1 the existence of z is proven. 

(3) One could adopt the point of view by Green, Rivlin and Naghdi along the 
path indicated by Marsden's and Hughes's theorem in [16] by postulating 
the expression of the first principle of thermodynamics (a point of view 
exploited in [3] with reference to isometric changes in observers). In this 
case, however, one is forced not only to postulate the existence of the energy 
but also to prescribe its functional dependence on the state variables, in 
contrast with the minimalist approach followed in Theorem 1. Such a point 
of view is however one of the manners useful to prove the covariance of the 
pointwise balance of actions (that is the invariance under the action of the 
entire group of diffeomorphisms on the ambient space and the action of G 
on A4). The other ways are given by the exploitation of Noether theorem 
and/or d'Alcmbert-Lagrange type principles in presence of viscous-type 
dissipation at macroscopic and/or substructural level (see [4] and [13] for the 
relevant results). A requirement of covariance allows one to eliminate the 
indetermination given by zi G KerA*: in this case z-i vanishes identically. 

When v represents only a generic property of the substructures not associated 
with their geometry in space, changes of frame in the ambient space and over A4 
can be considered disconnected. By imposing invariance of the external power 
with respect to isometric semi-classical changes in observers, one gets two distinct 
integral balances of moments: 



I ' {y-y )Abdx+ ( (y- y ) A Pn dH^ 1 = 

Jb Jdb 
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and 

/ A* (3 dx + [ A*Sn dH d - 1 = 0. 

Jb Jdb 

Theorem 1 can be rewritten. The sole difference in this case is that 13.71 splits in 
the two equations 

(3.10) skw(PF*) = 0, skw{A*z + (DA*)S) = 0. 

4. The energetic scenario and the existence of ground states 

4.1. A priori restrictions on constitutive structures. After describing the 
morphology of the generic material element and representing the actions along a 
motion, the local energetic scenario must be specified. 

At the macroscopic scale, since deformation is accounted for, each material ele- 
ment (considered as a whole) is assumed in energetic contact with the neighboring 
fellows. The consequent interactions are standard tensions. 

At the scale of the substructure, i.e. within each material element, some al- 
ternatives are possible (see additional remarks in [14 ]). They are classified under 
suggestion of the common path followed in statistical physics. 

(1) The generic material element is a closed system with respect to its substruc- 
ture: there is no migration of substructures out of the material clement, and 
the substructure itself does not interact energetically with the neighboring 
fellows. 

(2) The substructure of the generic material element is in energetic contact 
with the substructures of the neighboring elements. No migration occur. 

(3) The material element is an open system: both energetic contact and mi- 
gration of substructures are possible. 

Here the attention is focused on case 2 mainly. Remarks are added on case 1. 
Case 3 is not touched here for the sake of brevity (see [T2] for relevant developments). 

The procedure to establish a priori constitutive restrictions is the standard one, 
based on the Clausius-Duhem inequality which is written here in isothermal form 
as a mechanical dissipation inequality. It prescribes that 

(4.1) |* ( b ) _ v e b xt (y,i>) < 0, 

for any choice of the rate fields. \& (b) is the total free energy of b along (x, t) i — > 
(y (x, t) ,v(x,t)). The standard assumption is that \f (b) is absolutely continuous 
with respect to the volume measure so that there is a density ip such that 

*(b) = [ ip dx. 

Jb 

The constitutive dependence on the state variables must then be assigned not only 
for ip but also for the stress measures (namely P, z and S). Simple assumptions 
are as follows: 

ip = ip(F,v,N), P = P (F, v, N) , z = z(F,v,N), S = S (F, v, N) . 

If ip admits partial derivatives with respect to its entries, the arbitrariness of b and 
equation (|3.9[) imply the local dissipation inequality 

(4.2) {d F ip -P)-F+ (d v ip -z)-u+ (d N ip - S) ■ N < 0. 
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The possibility to choose arbitrarily the rate in (|4.2|) from any given state (F, u, N) 
implies the classical relations (see also [T]) 

(4.3) P = d F ip(F,v,N), z = d v iP{F,v,N), S = 8 N yj (F, v, N) . 

The mechanical dissipation inequality (|4.1| forbids the dependence of i/> on the 
rate of the fields involved. In fact, if -0 would depend on (let say) v, in the reduced 
version of the mechanical dissipation inequality a term of the type d^ip (F, v,N,v)- 
v would appear with no correspondence in the structure of the internal power 
VI (y, v) where no action developing power in v is presented. The arbitrariness 
of v would imply then d^tp = 0. In contrast, P, z and S may depend on the 
rates of the state variables when viscous-like effects occur at various scales. The 
dependence on the the rate of the state variables is compatible with the mechanical 
dissipation inequality 14.11 provided that one assumes the validity of an additive 
decomposition of P, z and S into conservative and dissipative parts. By indicating 
the triple (F, is, N) by one then presumes that 

P = P c (0 + P d M), z = z c (,) + z d (,,i), S = S C (?) +S d (?,<;). 

Such decompositions must be supplemented by the assumption that the conserva- 
tive components are determined by the free energy. The use of (14. ip implies once 
more the relations (|4.3[) for the conservative addenda and the reduced dissipation 
inequality 

(4.4) P d -F + z d -v + S d -N > 0. 

Consequently, P d , z d and S d are linear in F, i> and N. Additional assumptions on 
the structure of the dissipation can be made. 

(1) One may presume that strong dissipation conditions are satisfied a priori 
(that is independently of 14. ip . namely 

P d ■ F > 0, z d ■ v > 0, S d ■ N > 0. 

Then one may write 

P d = a P F, z d = a z p, S d = a s N, 

with ap, a z and as positive definite (scalar valued) state functions. 

(2) One could consider a strong condition for the macroscopic dissipation, 
namely 

P d ■ F > 

and a weaker dissipation condition for the substructure: 

z d ■ V + S d ■ N > 0. 

In this case, P d is equal to apF while z d and S d are linear functions of i> 
and N. 

(3) Other conditions can be presumed to hold. They may describe different 
viscous-like effects. Dissipative effects of plastic-like type can be accounted 
for. The standard plasticity theory and its strain gradient version fall within 
this scheme, when one identifies v with the plastic strain. In all cases of 
plastic-like behavior a flow condition in terms of the subdifferential of some 
admissible region in the state space must be involved (here the existence of a 
dissipation pseudo-potential is assumed) . Energetic solutions to the result- 
ing evolutionary problem can be obtained under appropriate hypotheses 
(relevant analytical tools can be found in [17], [6]). 
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One may ask what is the relation with standard internal variable schemes. In 
their initial formulation, such schemes have been proposed with the aim of describ- 
ing the removal from thermodynamical equilibrium (see [5] , [10] V In this (historical) 
sense internal variables are by definition not observable and play a parametric role 
at equilibrium, in contrast with the approach proposed here. The derivatives of 
the energy with respect to the internal variables and their derivatives are not con- 
sidered true interactions, rather they are thermodynamic affinities (see once more 
[17j . [5]). They do not appear in the expression of the external power. In contrast, 
I consider v as an observable quantity the variations of which contribute to the 
equilibrium by means of true interactions. This is the reason for which I call v 
morphological descriptor rather than internal variable. Connections are possible 
between the internal variable scheme and the multifield scheme that I discuss here. 
Assumptions should be necessary in order to avoid to render the comparison only 
formal. 

4.2. Ground states. Consider a complex body displaying a pure conservative be- 
havior. In this case one may identify the free energy with the elastic one. In absence 
of inertia, the energy of the whole body is then 

£{v, v ) : = / e(x,y{x),F(x),f(x),N(x)) dx, 

JB 

where the density e is the difference between the elastic energy e % and the potential 
of body forces ef + e|, the latter being decomposed in the part associated with 
standard gravitational forces (ef ) and the potential of possible external fields acting 
directly over the substructure (e|), namely e = e % {x,F,v,N) — {e\ (u) + e| {v)), with 
e l the elastic energy. 

A pair of fields (it, v) satisfying the variational principle 

min£ (y, v) 

is called ground state. Conditions for the existence of ground states follow consti- 
tutive assumptions on (i) the nature of the functional classes in which one places 
y and v, (ii) the 'structural' properties of e. 

Let y : B — > R d be a Sobolev map, namely an element of W 1,1 (B,R d ). Denote 
first by M (F) the d— vector in M. 2d collecting all the minors of F (i.e. of Dy). 
M (F) is then an element of A d (R d x R d ). 

It is possible to construct the d— current integration G v over the graph of y. 
Precisely, G v is the linear functional on smooth d— forms u> with compact support 
in B x R d defined by 

G y := / (u(x,y(x)),M(Dy(x))) dx. 

JB 

The boundary current associated with G y is indicated by dG v and defined by 
dG y (u) := G u {duj),uo E E> d - 1 ( J BxR d ) withD d - 1 (BxR d ) the space of -forms 
with compact support in B x M. d (details on the nature and the properties of Carte- 
sian currents can be found in [S]). 

The functional spaces in which the existence of minima is investigated must be 
specified. Their choice has constitutive nature. 

(1) The macroscopic deformation y is assumed to be a weak diffeomorphism 
(in symbols y G dif 1 ' 1 (B,R d )). In fact, y is considered a W 1 ' 1 (B,R d ) map 
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such that (i) \M(Dy)\ € L 1 (B), (ii) dG y = on V d ~ l (B x R d ), (in) 
det Dy (x) > for almost every (iv) for any / S C^°(B x R d ) 

/ / (a;)) det (x) dx < sup f(x,w) dw. 

Jb JR d xeB 

In particular, the subspace 

dif r ^(B,R d ) ■= {yedzf^(B,R d )\ \M (Dy)\ e V , 

for some r > 1 , is of special interest below. 
(2) It is assumed that (i) M. has Riemannian structure with (at least) C 1 — metric 
gMi and (i*) covariant derivatives are explicitly calculated by making use 
of the natural Levi-Civita connection. The C 1 — Riemannian structure im- 
plies that M. can be isometrically embedded in a linear space isomorphic 
to K M (for some M) by Nash theorem: it is then considered as a closed 
submanifold of R M . It is then assumed that v 6 W 1,s (Bq,M.), s > 1, with 

iy M (B, M) := {v £ W M (£>, M M ) | 1/ (a?) € M for a.e. x} . 

The energy functional £ is then extended to 

W Pl , :=(ft/ r ' 1 C8,I d ) X W M (B..M). 

Assumptions on the structural properties of the energy density e must also be 
specified, e cannot be convex in F for the standard objectivity argument but it can 
be convex in N. 

• e is assumed to be poly convex in F and convex in N. There exists a Borel 
function Pe : B x R d x M x A d (R d x K d ) x M Mxd -> K+, with values 
Pe (x, y, f , £, iV), which is (i) 1. s. c. in (y, v, £, AT) for a.e. ifB, (ii) con- 
vex in (£, A) for any (x, y, v), (Hi) and also such that Pe (x, y, v, M (F) ,N) = 
e (x, y, v, F,N) for any (x, y, v, F,N) with det F > 0. 

• By assumption e satisfies the growth condition 

e (x, y, v, F,N) > C\ (\M (F)\ r + \N\ S ) + (det F) 

for any (x, y, v, F,N) with det F > 0, r, s > 1 and C\ > constants, and 
i? : (0, +oo ) — > M + a convex function such that (t) — > +oo as t — > + . 

Theorem 2. (|15j ) TTie functional £ achieves the minimum value in the classes 

W d s := {(y, v) e Wr, s \y = y a on dB y , v=v on dB v } 

and 

W r % := {(y, v) e W r , s | SGj, = dG Vo on V 2 (R d x M d ), i/ = i/ on dB v } . 

In the theorem above dB y and c?i?^ are the portions of the boundary dB where 
boundary data are assigned in terms of y and v respectively. Details and comments 
on the physical consequences of the assumptions above can be found in [15] . 
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5. Notes and complements 

The first variation of the energy functional £ (y, v) along C 1 minimizers allows 
one to obtain the balances of standard and substructural actions (|3.6p and (|3.8[) . 
The condition f|3 . T|) is a consequence of a requirement of objectivity for the elastic 
energy e l . Remarks leading to (13.10|) also apply. In addition, horizontal varia- 
tions can be made by altering the reference place by means of the diffcomorphism 
$ £ (x) := x+£(f> (x) , (j) € Cq (£>, R 3 ) , e a real parameter. $ e (x) leaves unchanged the 
boundary dB for e sufficiently small. In fact, horizontal variations can be considered 
as a sort of relabeling of the reference place. One then defines y e (x) :— y (^J 1 (a;)) 
and u £ (x) := v (^ 1 (x)), and obtains a mapping e i — > £ (y e , v e ). In case of appro- 
priate smoothness, differentiation with respect to s gives rise to the configurational 
balance 

(5.1) DivP + d x e = Q. 

where P = e l I - F*P - N*S G Aut (R d ) is the extended Hamilton-Eshelby tensor 
in the mechanics of complex bodies (see [H]). In the smooth case (|5.ip is essentially 
the pull back in B of the balance of standard forces (|3.6p . 

Different is the stage when the variation is calculated on non-smooth minimizers, 
namely on Wf s where it is sure by Theorem 2 that ground states exist. The 
main difficulty is that Sobolev maps may not admit tangential derivatives so that 
the balance of standard forces in terms of first Piola-Kirchhoff stress (|3 .6[) cannot 
be derived. One may compute horizontal variations, variations over M. and the 
variation of the actual shape B a . I summarize here the relevant results (see for 
details ft5]). 

(1) Under appropriate growth conditions for the polyconvex extension of the 
energy and its derivatives with respect to x, M (F) and N, one proves that 
(i) F*P and N*S belong to L 1 (B) and (ii) the balance (j5.ip holds in terms 
of distributions. 

(2) An additional assumption on the growth of \d y Pe\ implies that (i) a := 
(detF) -1 PF* e L\ oc {y (B ) ,R 3 ® R 3 ), with x\ — > y {x) the Lusin repre- 
sentative of y and (ii) the balance of standard forces hold in distributional 
sense in terms of the Cauchy stress a. 

(3) Variations over AA, obtained by means of smooth curves e — > (p s <E Aut (Ai), 
(p G C 1 (M), and an additional assumption on the growth of \d v Pe\, allow 
one to show that S G L 1 (B ,M. 3 * ®T*M) and (ii) the balance of sub- 
structural actions (|3.8|) holds in the sense of distributions. Such a result 
excludes the interpretation (suggested by some author) of the balance of 
substructural actions as a sort of balance of configurational actions, unless 
the word configurational is used in a sense a little bit far from the current 
one. 

The remarks in this section contribute to the current debate on the nature of 
configurational actions. 
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